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Arbitrary waves incident on a solid embedded nanoparticle are studied. The acoustic vibrational frequencies 
are shown to correspond to the poles of the scattering cross section in the complex frequency plane. The 
location of the poles is unchanged even if the incident wave is nonplanar. A second approach approximating 
the infinite matrix as a very large shell surrounding the nanoparticle provides an alternate way of predicting 
the mode frequencies. The wave function of the vibration is also provided. 



1 Introduction In the last two decades, a significant number of works have been devoted to the study 
of acoustic vibrations of spherical nanoparticles. This has aided the understanding of the role of vibra- 
tions in the performance of some optical devices (for example, in electronic dephasing due to emission of 
phonons). These vibrations can be experimentally investigated with different techniques including Raman 
scattering 1 1 1, photoluminescence|2|, and femtosecond pump-probe experiments! 3 1. Usually, continuum 
elastic models have been used to describe these vibrations. The original 1882 work by Lamb[4| found the 
vibrational frequencies of a free continuous isotropic elastic sphere. However in real life, nanoparticles are 
seldom free but rather embedded in another medium such as a glass matrix or liquid. 

The problem of eigenmodes of vibrations of an embedded sphere has not been extensively studied. 
However, a lot of work has been devoted to the apparently unrelated problem of scattering of acoustic 
waves by a sphere|5 1. That is where Dubrovskiy et al. paper|6| has its starting point to calculate complex 
natural frequencies of vibrations of an embedded sphere . Although this method, referred to here as the 
complex frequency model (CFM), is a very powerful tool, one might wonder about the significance of the 
complex frequencies and the relevance of discrete natural vibrations in a problem where one expects a 
continuum of vibrational state. The goals of this paper are to explain the basics of this model and to show 
its relevance by comparing it with the core shell model (CSM)|7| which is a different approach using real 
valued frequencies. 



2 Scattering of elastic waves and CFM The usual basis of unit vectors in spherical coordinates is (e r , 
eg, e^). However, it is more convenient to express results using a different set of vectors. In this paper, we 
will use basic vectors (e\, e2, £3) defined as follows: 

ei = Y e m e r (1) 

1 dY ( m A drr . 

sin6> oq> oO 

1 drr A drr „ 

sin 6* dm 30 



* Corresponding author: e-mail: lucien.saviot@U-bOurgogne.fr 



Copyright line will be provided by the publisher 



2 



L. Saviot 1 and D. B. Murray 2 : Acoustic vibrations of an embedded nanoparticle 



where Y/ 11 represents the spherical harmonics. 

Vibrational modes of a sphere, whether free or embedded, may be classified as "torsional" (TOR) or 
"spheroidal" (SPH). Torsional modes are those for which the displacement field has zero divergence. All 
other modes are called spheroidal. Torsional modes have non-zero displacement along e 3 only. Spheroidal 
modes have non-zero displacement along e\ and e 2 only. For £ = 0, we have £2 = £3 = 0. 



2.1 Basic equations In order to clarify the meaning of the complex-valued eigenfrequencies of an 
embedded elastic sphere, the problem we want to investigate is the scattering of an incident elastic wave 
in a medium (matrix) by a sphere made of a different material than the matrix. Three different terms have 
to be considered: the incident wave, the motion inside the sphere and the scattered wave. All these waves 
have to satisfy the following wave equation (assuming an exp(— iuit) time dependence) inside the medium 
where they exist. 

(A+2 ) u)V(V-m) -/iVx(Vxu) + puj 2 u = (4) 

u is the displacement, A and a the Lame coefficients and p the mass density. 
Thus, they can be decomposed into the the following form: 

u = V$ + Vx(rx + Vx(r0)) (5) 
where the three potential fields are 

$ = A im z e {k L r)Y e m (6,<t>) (6) 

l=0,m 

X = Yl B emZi (k T r)Y e m (6,cb) (7) 

l=l,m 

4> = J2 Ce m ze(k T r)Y e m (0,cf>) (8) 

l=l,m 

zg is the relevant spherical Bessel function of order I, (first or second kind or Hankel) k^ and hr are 
longitudinal and transverse wavevectors respectively. 



2.2 Solutions Using the notations introduced by Eringen et al. | 8 1, the displacement and the associated 
surface traction (surface force per unit area) T = a rr e r + a r g eg + oy^ of the (£,m) component can be 
written as follows: 



potential 
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(9) 



The z superscript has been added to recall the nature of the spherical Bessel function to use. Inside the 
sphere, z is the spherical Bessel function of the first kind j. Outside the sphere, in the general case, two 
terms have to be considered where z is j and n (spherical Bessel function of the second kind). For the 
case of the outgoing scattered wave, we'll use the spherical Hankel function of the first kind h/p{x) = 
ji(x) + int(x). 
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2.3 Boundary conditions The boundary conditions are the usual continuities of displacement and sur- 
face traction at the sphere-matrix interface. Because the boundary conditions should be satisfied for r = R 
(sphere radius) at all 6 and <f>, each equation of continuity breaks up into a linear inhomogeneous system 
of six equations for each (£, m) the incident wave is decomposed onto. The unknowns are the amplitudes 
of the vibration inside the sphere and of the scattered wave. The constants are the amplitudes of the de- 
composition of the incident wave on (£, to). Therefore the unknowns can be expressed as the ratio of two 
6x6 determinants. The determinant equation in the denominator is the one of the homogeneous system, 
i. e. fori > and |m| < t 
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(10) 



where £ = UlR, r\ = krR, H P and \i m are the particle and matrix mass densities and the rows represent the 
continuity of the displacement along u e\, u-e 2 , T-e±, T-e 2 , u e^ and T-e% from top to bottom. The columns 
represent the coefficients for $ p , ip p , $ m , tp m , \p, and \ m where the p and to subscripts are for waves 
inside the particle and the outgoing wave inside the matrix respectively. 

Since the matrix is block diagonal this determinant is easily decomposed as the product of A; (4 x 4 
determinant for SPH modes) and ^(2x2 determinant for TOR modes). 

For £=0, both e 2 and C3 equal and the determinant is: 



D = A = 



U((£ = 0,0 Uf\£ = 0,0 



(11) 



These determinants do not depend on the azimuthal quantum number to. They are also independent of 
the nature of the incident wave. This makes the results equally applicable to a situation where the incident 
wave is a spherical wave emanating from a source close by to the scattering sphere. This would be the case 
in a phononic crystal, where elastic waves scatter sequentially from a periodic array of scatterers. 



2.4 CFM Up to now, only real valued uj have been considered. The amplitudes of the wave inside the 
sphere and of the outgoing wave have been shown to be the ratio of two determinants with the same de- 
terminant in the denominator (the determinant of the homogeneous linear system). These amplitudes can 
therefore become quite large when this determinant is close to zero. The physical system corresponding 
to the system of homogeneous equations is the one where there is no incident wave. In this case, the only 
solution is obtained when all the amplitudes are zero. Therefore Di does not vanish for real values of to. 
However, it can reach local minima. In order to predict the position of these minima, we are going to locate 
the complex roots of the matrix determinant D(. This is the heart of the complex frequency model intro- 
duced by Dubrovskiy and Morochnik|6|. Because of the cxp(— icut) dependence of the displacement, the 
real parts of the CFM roots correspond to the position of the resonance and the imaginary parts correspond 
to its width and its damping. Of course, this interpretation is just an approximation and we can't expect it 
to precisely reproduce the position and width of the embedded sphere resonances in every case. As seen 
before, A; does not depend on the nature of the incident wave. Therefore its roots can be associated with 
the embedded sphere natural vibrations which we prefer to term pseudo-modes of vibrations. 

Two kinds of CFM roots exist. Some reflect primarily natural vibrations of the sphere. Some others are 
much more damped and reflect primarily the vibrations of the matrix around the sphere. The latter (later 
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referred to as "matrix modes") appear as a background in the scattering calculations. Sphere vibrations 
are narrower peaks above this background. Of course, this distinction becomes harder to make when the 
sphere and the matrix have similar characteristics (mass density and sound velocities). Therefore, there are 
several limitations inherent to the CFM model: 

1 . it is an approximation (because knowing the complex roots of A/ is not enough to know its variation 
for real arguments) 

2. what's the relative importance of sphere and matrix modes? 

3. this approach does not provide valid wave functions which are required for Raman scattering spectra 
calculations for example 

3 CSM To overcome the limitations of the CFM model, it is necessary to use a different approach. 
Instead of starting from a formulation originating from an elastic wave scattering problem, it seems more 
interesting to focus on the specificities of optical properties of nanoparticles. In the Raman scattering prob- 
lem, the vibration eigenmodes and eigenenergies of the system are important. Then, for a given eigenmode, 
the Raman intensity is related to the displacement inside the particle (in most usual cases, Raman scattering 
from the matrix is negligible in this energy range). In a recent paper|7 1, we modeled the matrix embedded 
nanoparticle by a core-shell system where the "core" is the nanoparticle and the "shell" is the matrix. The 
shell's external radius is chosen several orders of magnitude bigger than the nanoparticle radius. The mean 
squared displacement within the nanoparticle interior ((u 2 ) p ) was monitored by integrating the squared 
displacement field over the interior of the core as a function of the eigenmode energy. This uses displace- 
ment fields that were normalized over the entire core and matrix. As was shown, very good agreement is 
obtained between the frequency and width of particle modes calculated with CFM and peaks in the plots 
of (u 2 ) p as a function of the mode frequency. The validity of the CFM model is challenged when the 
nanoparticle and the matrix are made of materials having similar mass densities and sound velocities, for 
example for a silicon nanoparticle inside silica. In such cases, the CFM particle modes are still in good 
agreement with the peaks in the (u 2 ) p plots, but a significant background is seen below these peaks. 

4 Conclusion Depending on the precision required for a given task, there are several ways one can 
approach calculations of natural vibrations of a free or embedded elastic sphere. If the sphere is free, then 
the Lamb model suffices. When the sphere is not free but only the frequency of the natural vibrations 
matters, the Lamb model (free sphere model (FSM)) or the bound sphere model|7| (BSM) will suffice in 
most cases as a workable approximation. When the width of the resonance matters the CFM is an accurate 
tool. CFM also fills in a few special cases not addressed by FSM, such as the existence of matrix modes. 
Finally, for situations such as the calculation of the absolute intensity or shape of low-frequency Raman 
spectra, models such as the CSM are necessary. 
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